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On  the  Negative  Binomial  Convergence 
in  a Class  of  m-Dimensional  Simple  Epidemics 

by 

H.  Lacayo  and  Naftali  A.  Langberg 
ABSTRACT 


We  consider  a population  which  is  exposed  to  m infections,  and 
consist  initially  of  N susceptibles.  At  each  point  in  time  at  most  one 
susceptible  becomes  infective,  and  only  from  one  cause.  This  m-dimensional 
simple  epidemic  is  a stochastic  process,  (X  (t),  ....  X..  (t)),  with  conponents 

counting  the  number  of  infectives  from  the  respective  causes  at  time  t. 

We  show  that  if  the  transition  rates  of  cause  1 through  m at  tine  t are 


, m,  a. 


a >0 
m 


give*  by  YjCHI  - ‘ I a,  (t)  - X,  (O)J-I  i - 1, 

1=1 

and  if  Lim  X^,  . (0)  = b.  e(l,  2,  then  (X.  . (t) , ....  XM  (t))  converges  as 

N 4 • to  a random  vector  with  independent  negative  binomial  components. 


Key  Words:  Exponential  distributions,  convergence  in  law,  stochastic 

processes,  m-dimensional  simple  epidemics,  negative  binomial 
distributions. 
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1.  Introduction  and  Summary 


Introduction. 

We  consider  a population  which  is  exposed  to  m infections  and 
consists  initially  of  N susceptibles  and  bM  . infectives  from  cause 
i,  i ■ 1,  ...»  m.  At  each  point  in  time  at  most  one  susceptible  becomes 
an  infective,  and  only  from  one  cause.  Once  an  individual  enters  the 
infective  state  he  remains  there  and  cannot  be  a carrier  of  any  other 
infection.  This  m-dimensional  simple  epidemic  is  a stochastic  process, 

(*N  1^5*  •••»  m (•£)),  v/ith  components  counting  the  number  of  infectives 
from  the  respective  causes  at  time  t.  We  show  that  if  the  transition 
rates  of  cause  1 through  m at  time  t are  given  by: 


■iVitt),;i  -w  ’Vi”- *-  *■ 

1=1 


i m,  a^i  *..,  ^ 0 , 


(1.1) 


and  if 


Lira  bM.«b.  e(l,  2,  ...},  i=l,  ...,m  (1.2) 

xt ■-  ^ » 1 1 

IT**0 


then  for  every  positive  real  numbers  t and  6,  and  m nonnegative  integers 
kl»  km' 

m 

Lim  POL,  .(t)  a.  +b.,  i«  1,  ....  m)  * n P{X.(t)  a.  ♦ b.)  (1.3) 

N~  1 1 i-1  1 1 1 


and 


E(Xi(t))8,  i - 1,  ....  m. 


where  Xj (t) , , 


(1.4) 


X (t)  are  negative  binomial  random  variables  with 
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respective  parameters  b.  and  e ^ , i - 1,  ....  ».  The  only  previous  work 
in  the  area  of  m-dimensional  simple  epidemics  is  the  one  by  Billard,  Lacayo 
and  Langberg  (1978),  who  proved  (1.3)  and  (1.4)  with  the  additional 
assumption  that  «1  ■ a2  » am.  In  that  case  the  interinfection 
random  times  are  independent  of  the  infection  causes.  Generally,  as 
will  be  pointed  out  in  Section  2 this  property  does  not  hold.  The  proof  of 
(1.3)  in  the  cited  reference  depends  on  the  special  structure  of  the 
interinfection  times.  To  obtain  (1.3)  and  (1.4)  we  have  approached  the 
problem  from  an  alternate  viewpoint. 

Summary. 

In  Section  2 we  present  a rigorous  definition  of  an  m-dimensional 
simple  epidemic,  and  describe  the  ones  that  are  the  subject  of  our 
analysis.  Statements  (1.3)  and  (1.4)  are  proven  in  Section  3.  In  this 
section  we  present  an  m-dimensional  simple  epidemic  in  a population  con- 
sisting initially  of  infinitely  many  susceptibles.  This  theoretical 
simple  epidemic  is  instrunental  in  the  proofs  of  (1.3)  and  (1.4). 
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2.  The  m-dimensional  Simple  Epidemic 
An  m-dimensional  simple  epidemic  in  a population  consisting  initially 
of  N susceptibles  and  t>N  ^ infectives  from  cause  i,  i * 1,  . ..,  m can  be 
described  by  N bivariate  random  vectors  (T^,  CN>1),  ...,  (TN>N,  CN  N). 

TN  1’  Tn  N’  are  the  random  interinfection  times,  and  CN  , , Cu  M 

are  discrete  random  variables  with  values  in  {1,  m}  designating  the 

k 

respective  infection  causes.  Let  SN  ..  be  equal  to  £ T , k * 1 , . . . , N. 

* j*l  * 

The  finite  dimensional  joint  distributions  cf  (X.,  , (tl , . . . , x . ft))  are 

N,1 v - * * N,m 

determined  by  (TN  j,  j) , ....  (TN>N»  ?N  N) » through  the  following  set 
equalities: 

Yi(t)  * Yi  * k)  * (SN,r  5 * ‘ SN.r.l-  Z.^Yq  " «>  " k>.  * » « C2-» 

T— K q«l  n 

k - 0,  N,  and  i - 1,  ....  m.  (SN  Q = 0,  SN  N+1  = -) . 

We  assume  throughout  that  (T^j,  CN>1),  ...  (Tj.  N,  ?N  N) , are  determined 
by  (2.2)  and  (2.3). 


PUN,j  * ly  3 * i»  •••.»> 


[ n a *-  ni_1(b  ♦ j)X  n ( l a.(b  ♦Tk*  * i))]"1, 

i-1  1 j-0  N'1  k«l  i-1  1 N>1  ml  3 

N 

where  r • J 1(1  • i),  i « 1,  ....  m. 

rl  • N 

z ® 

P{TN,k^N,l*  5N,k-l  * V k * *»  * e j-1 

where  tJ(  . . . , tR  are  positive  real  numbers,  and  p,  is  equal  to 

n - j ♦ i 5 , i;1  5 

N I «i (bi  ♦ [ Itta  ■ i)),  for  J ■ 1,  ....  N. 

i-1  1 1 c^l  q 


N m 


k-1 


(2.2) 


(2.3) 


In  the  particular  case  discussed  in  Billard,  Lacayo  and  Langberg  (1978), 
condition  (2.3)  reduces  to  independent  exponential  interinfection  times 


of  Cu  ,,  . . . , In  addition  equation  (2.2)  simplifies  to 

m N4i-1  N>N  N m , 

( n n (b„  . + j)][  I ( [ b . + k - 1)1.  Finally  we  note  that  if 
i«l  j»0  ’ k=l  i= 1 

(2.2)  and  (2.3)  hold,  then  the  m-dimensional  simple  epidemic  satisfies 
(1.1).  Conversly,  if  (1.1)  is  satisfied,  then  equations  (2.2)  and  (2.3) 
follow  by  the  inherited  Markovian  structure  of  the  m-dimensional  simple 
epidemic. 
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3.  Main  Results 

Let  (Xjft),  ....  Xm(t))  be  an  m-dimensional  simple  epidemic  in  a 

population  consisting  initially  of  b ^ infectives  from  cause  i, 

i * 1,  2,  . ..,  m,  and  infinitely  many  susceptibles.  This  epidemic  can  be 

described  by  a sequence  of  bivariate  random  vectors,  (T^ , £^),  (T2,  £2) . 

....  2,  ...  . Tj,  Tj,  ...»  are  the  interinfection  times,  and 

£1*  €?,  •••»  are  random  variables  with  values  in  {1,  ...,  m},  which 

k 

designate  the  respective  infection  causes.  Let  S,  be  equal  to  J T., 

3 

k * 1,  2,  ...  . The  finite  dimensional  joint  distributions  of  (Xj(t), 

Xm(t))  are  determined  by  (T^,  5^) , k « 1,  2,  ....  through  the  following 
set  equalities: 

00  X* 

(X  (t)  = b.  ♦ k)  - U(SrSt<  s l I(c  = i)  = k)  (3.1) 

r*=k  q*l  4 

for  every  positive  real  number  t,  positive  integer  k,  and  i * 1,  ....  m. 

We  assume  throughout  that  the  sequence  (T, , ) , (?2,  £2) , ...,  is  deter- 

mined by  (3.2)  and  (3.3). 


3 ^ •••»  n} 


m r,-l 


[ n n*  *(b  ♦ j)JC  n ( f a.(b.  * Jfi(*  = djt 

i-l  j-0  1 j«l  i-1  11  q=l  *1 

n 

where  r.  ■ £ I(i  • i) , i = 1 , . . . , m,  and  n * 1,  2,  . 


(3.2) 


n 

- I 


WTfcHi.  •••»  > tk,  k ■ 1,  ....  n)  = e j-1 


i 3 


(3.3) 


r J;1 

where  tj,  ...»  tn  are  positive  real  numbers,  p.  ■ Z c^fb.  * Z ■ 1)) 

■*  i»i  q«l  ^ 

j * 1,  . . . 1 n , and  n ® 1,  2,  ...  . 
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We  will  show  that  for  every  two  positive  real  numbers  t and  B 


xm  „,(*))  converges  in  law  to  (X.  (t),  X (t))  as 

in, m i m 


N -*•  and  that  Lim  E(XM  . (t))B  = E(X. (t))B  for  i = 1,  ....  m.  First 
we  introduce  two  lemmas. 


Lemma  3.1.  Let  Uj,  U2,  ...,  be  an  i.i.d.  sequence  of  exponential  random 

variables  with  mean  1.  Further  let  (r.  t.  , . ,)  be  the 

D+r-i  1 D+r-l ,D+r-i 

order  statistic  of  a sample  of  size  b + r - 1 from  U. . Then  £ (b  ♦ j - l)"*u, 

j=l  3 

and  Tjj+r_j  r are  identically  distributed. 

Proof.  It  suffices  to  note  that  the  spacings;  t . , . , - t 

K 6 ’ r+b-l.r-j+l  r+j-l,r-j' 

j =»  1,  ...»  r,  0 = 0)  are  independent  exponential  random  variables 

with  means  equal  respectively  to  (b  ♦ j - 1)"*,  j = 1,  ...,  r. 

m 

Lemma  3.2.  Let  t and  B be  positive  real  numbers;  let  c = Sup  ][  b . and 

N i=l  N,i 

a » max  a. . If  (1.2)  holds  then;  (i)  P(SW  s t)  5 (1  • e"“ct)r  for 
isism  1 „ ' 

r a 1,  .... 


N,  and  (ii)  Sup  E(XN  . (t))p  < »,  for  i = 1,  ...,  m. 

•1  IN*1 


Proof.  We  note  that  (ii)  is  a corollary  of  (i).  To  prove  (i)  let 

r , 

Uj,  U2,  ....  be  as  in  Lemma  3.1.  Since  P{SN  ^ s t)  £ P{  £ j Ih  s oct), 
(i)  follows  from  Lemma  3.1. 


Theorem  3.3.  Let  8 and  t be  positive  real  numbers.  If  (1.2)  holds  then; 

(i)  •••»  X^m(t))  converges  in  law  as  N + - to  (X^t),  X^t)) 

and  (ii)  Lim  E(XNi(t))6  - E(X.(t))B  for  i a 1,  ....  n. 
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Proof.  We  note  that  (ii)  is  a corollary  of  (i)  and  Lemna  3.2  (ii). 

To  prove  (i)  let  k^,  ...»  be  m nonnegative  integers,  A be  the  set 

m 

{(rj,  ....  j r^^  a kj,  i * 1,  ....  ra}  and  let  R =>  £ ri#  Since  for  N 

sufficiently  large,  P{X  (t)  a k.  + b. , i = 1.  . . . , m}  = 

fr  1 

Z^^N,R  5 * < ^N,R+1  ’ ^ ~ Ti»  - = l»  •••»  m)I(R  - H) } , (i)  follows 

from  Lemma  (3.2)  (i) , the  dominated  convergence  theorem,  and  some  simple 

calculations. 


For  reference  purposes  the  following  two  well  known  results  are  presented. 

Lemma  5.4.  Let  U^,  Uj,  . . . , be  as  in  Lemma  3.1,  and  let  p^,  v»2 » •••»  he 
a sequence  of  positive  real  numbers.  Further  we  denote  by  f^  the  density 


function  of  £ p. *U. , k = 1,  2,  . . . . Then,  (i)  p-1.  f . (t)  * 

r j=l  \*13  r T+1r  T+A  r , 

P{£  p?V  st  < l pT  U },  and  (ii)  f ft)  = ( IT  u.)  I e'^j1  n (p.  - p.)'\ 
j=l  3 3 j=l  3 3 T j-1  3 j=l  i~l  3 

provided  p^,  ...,  pr  are  distinct.  x*3 


-1 


We  are  ready  to  state  and  prove  the  main  result  of  the  paper. 

Theorem  3.5.  Let  t and  B be  positive  real  numbers,  and  k, , k be 

l m 

nonnegative  integers.  Further  let  A be  the  set  (r^,  ...,  r^jr^  z k^, 
i ■ 1,  ...»  m.}.  If  conditions  (1.1)  and  (1.2)  are  satisfied,  then  the 
following  hold. 

m fk  x w _ 1 1 


m 


(i)  Lim  POL,  (t)  a.+b,i»l(  ....  m)  - J n 
N+«®  * A i=l 


b.  + r.  - 1 

l l 


bi  ' 1 


1 


e’ait)  \ 


, B 


(ii)  Lim  E(X..  .(t))  is  the  8th  moment  of  a negative  binomial  random  variable 

Um.  


-Oit 


with  parameters  b^  and  e 1 , i =»  1 


m. 


i 
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Proof.  By  Theorem  3.3  it  suffices  to  show  that  X (t) , ....  x (t)  are  inde- 

■i  m 

pendent,  and  have  the  respective  negative  binomial  marginal  distributions. 

Let  kj,  • ...  be  nonnegative  integers,  with  sum  equal  to  k.  Further,  let 
B denote  the  set 


CUj,  ...»  *k)Ur  ....  e {1,  ...,  m},  l I(£  - i)  - k.,  i ■ 1,  ....  ml. 

q=l  M 

We  note  that  P(X.(t)  = k.  4.  b.,  i = i,  . . . , m>  equals 

l P{Sk  5 t < Sk+1|?q  = *q,  q = 1,  ...,  k}  P{^  = £q,  q = 1,  ...,  k}. 

From  Lemma  3.4  and  equation  (3.2)  it  follows  that  P{X.(t)  = k.  ♦ b.,  i « 1,  ...,  n} 
equals 


m 

n 

i*l 


b.  + k.  - l3  -a.b.t  k. 


1 1 

b.  - 1 
1 


1 1 « H • 


k+1  -e.t  k+1 


i ki'  l I e ' n where 

B j=l  i=l  1 3 


j-1  m 

e-i  = I I Bjl(l  5 i),  j * 1,  .... 
J q=l  i=l  1 q 


k + l, 


By  reusing  Lemma  3.4  we  obtain  that  P{X.(t)  = k.  + b. , i * 1 , . . . , k}  equals 


in 

n 

i=l 


bi  * kt  - 1 
bi  - 1 


-a.b.t  k. 

. l 1 1,  , 

e « i*  * 


Ki  r Kt1  -1  k+1  , 

ai  ki!  Z p<  Z 04 \ st)  n e?1 

B j=l  3 3 j*2  3 


To  complete  the  proof  it  suffices  to  show  that 

m ki  r kt3  -1  k+l  1 m -a.t  k. 

n a kj!  I p{  l 8 .u.  st}  n e"1  = n (1  . e 1 ) i. 

1-1  B j-2  3 3 j=2  3 i=l 

Let  F denote  the  left  side  of  equation  (3.4).  Further  let 

It 

®i  * *’2'  4q^lA2’  *v  e (*»  •••*  j I(£  - r) 


q«2 


IV 

kr»  r • I.  ...»  10,  r * i,  £ l(*  - i)  . v - 

q«2  q * 


(3.4) 


13. 


m j-1 

I aT  1 I(*  ■ r),  i * 1,  . ..,  m,  j = 3,  . ...  k + 1.  Consequently  Equation 

r=l  q*2  ^ 

r*i 

(3.4)  is  established  by  an  induction  argument  on  k applied  to  each  of  the 
k+2  k+1 

sums  l P{  l e.  .U  s t)  n 0T1  , i = 1,  ....  m. 

B.  ■»  j-3  3,1 

Corollary  3.6.  Let  X^(t)  the  total  number  of  infectives  at  time  t.  If 
(1.1)  and  (1.2)  are  satisfied,  then  the  following  hold. 

(i)  X^(t)  converges  in  law  as  N + « to  a sum  of  m independent  negative 

-a.t 

binomial  random  variables  with  respective  parameters  b.,  e 1 , i ■ 1,  ...,  m,  and 

g 

(ii)  Lim  E(X  (t ) ) equals  to  the  8th  moment  of  the  limiting  random  variable 
N-**>  N 

given  in  (i) . 


In  t’<e  particular  case  discussed  in  Billard,  Lacayo  and  Langberg  (1978), 

XN(t)  converges  in  law  to  a negative  binomial  random  variable  with  parameters 
m _ 

J b.  and  e a . This  fact  played  a major  role  in  the  proof  of  (1.3)  for 
i-1  1 

the  symmetric  case  as  presented  in  the  cited  reference.  Finally  we  note 
that  it  is  of  interest  to  investigate  the  asymptotic  behavior  of  m-dimensional 
simple  epidemics,  that  do  not  have  the  Markovian  structure. 
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20.  ABSTRACT 

We  consider  a population  which  is  exposed  to  m infections,  and  consist  initially  of  N 

susceptibles.  At  each  point  in  time  at  most  one  susceptible  becomes  infective,  and  only  from 

one  cause.  This  m-dimensional  simple  epidemic  is  a stochastic  process,  (X..  . (t) , ...,  X.,  (t), 

N,  l N,m 

with  components  counting  the  number  of  infectives  from  the  respective  causes  at  time  t. 

We  show  that  if  the  transition  rates  of  cause  1 through  m at  time  t are  given  by 
l m 

«i*N,iWC1  - l (Xn  *(t)  - XN,i(0))l,i  ■ 1,  ...,  m,  Oj,  ....  am  > 0 and  if  Lint  XN  4(0)  • 
i-1  * ' N-*00  ’ 

bj  e (1,  2,  ...},  then  (XN  j^t),  ...»  X^  ffl(t))  converges  as  N -*■  » to  a random  vector  with 

independent  negative  binomial  components. 


